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CZ3 

' We consider a generalized class of Keller-Segel models describing the chemo- 

taxis of biological populations (bacteria, amoebae, endothelial cells, social in- 
sects,...). We show the analogy with nonlinear mean field Fokker- Planck equa- 
tions and generalized thermodynamics. As an illustration, we introduce a new 
model of chemotaxis incorporating both effects of anomalous diffusion and 
exclusion principle (volume filling). We also discuss the analogy between bi- 
ological populations described by the Kcller-Segcl model and self-gravitating 
Brownian particles described by the Smoluchowski-Poisson system. 
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QO , 1. Introduction 

Os- 

-^-j. ■ The name chemotaxis refers to the motion of organisms induced by chemical 

signals [1]. In some cases, the biological organisms (bacteria, amoebae, en- 
■ dothelial cells, ants...) secrete a substance (pheromone, smell, food, ...) that 

0© , has an attractive effect on the organisms themselves. Therefore, in addition 

to their diffusive motion, they move preferentially along the gradient of con- 
centration of the chemical they secrete (chemotactic flux) . When attraction 
prevails over diffusion, the chemotaxis can trigger a self-accelerating process 
5^ ' until a point at which aggregation takes place. This is the case for the slime 

mold Dictyostelium discoideum and for the bacteria Escherichia coli. This 
is referred to as chemotactic collapse. A model of slime mold aggregation 
has been introduced by Patlak [2] and Keller & Segel [3] in the form of two 
coupled differential equations. The first equation is a drift-diffusion equa- 
tion describing the evolution of the concentration of bacteria and the second 
equation is a diffusion equation with terms of source and degradation de- 
scribing the evolution of the concentration of the chemical. In the simplest 
model, the diffusion coefficient D and the mobility x of the bacteria are con- 
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stant. This forms the standard Kcllcr-Scgcl model. However, the original 
Keller-Segel model allows these coefficients to depend on the concentration 
of the bacteria and of the chemical. If we assume that they only depend on 
the concentration of the bacteria, the general Kcllcr-Scgel model becomes 
similar to a nonlinear mean field Fokkcr-Planck equation. Nonlinear Fokker- 
Planck (NFP) equations have been introduced in a very different context, 
in relation with a notion of generalized thermodynamics [4]. As far as we 
know, the connection between the general Kellcr-Scgcl model and nonlinear 
mean field Fokkcr-Planck equations has been first mentioned in Chavanis [5] 
and developed in subsequent papers (see [6] and references therein). This 
analogy makes possible to interpret results obtained in chemotaxis in terms 
of a generalized thermodynamics. At the same time, chemotaxis becomes 
an example of great physical importance for which a notion of (effective) 
generalized thermodynamics is justified. 

The standard Keller-Segel (KS) model has been extensively studied in 
the mathematical literature (see the review by Horstmann [7]). It was 
found early that, above a critical mass, the distribution of bacteria be- 
comes unstable and collapses. This chemotactic collapse leads ultimately to 
the formation of Dirac peaks [8-20]. Recently, it was shown by Chavanis, 
Rosier & Sire [21] that, when the equation for the evolution of the con- 
centration is approximated by a Poisson equation [10,14,17], the standard 
Keller-Segel (KS) model is isomorphic to the Smoluchowski-Poisson (SP) 
system describing self-gravitating Brownian particles. The chemotactic col- 
lapse of biological populations above a critical mass is equivalent to the 
gravitational collapse of self-gravitating Brownian particles below a critical 
temperature [22]. Assuming that the evolution is spherically symmetric, 
Chavanis & Sire [21-28] were able to describe all the phases of the collapse 
(pre-collapse and post-collapse) analytically in d dimensions, including the 
critical dimension d = 2. 

Recently, some authors have considered generalizations of the stan- 
dard Keller-Segel (KS) model. Two main classes of generalized Keller-Segel 
(GKS) models of chemotaxis have been introduced: 

(i) Models with filling factor: Hillen & Painter [29,30] considered a 
model with a normal diffusion and a density-dependent mobility x(p) = 
x(l — p/ctq) vanishing above a maximum density <jq. The same model was 
introduced independently by Chavanis [5,31] in relation with an "exclu- 
sion principle" connected to the Fermi-Dirac entropy in physical space. 
In these models, the density of bacteria remains always bounded by the 
maximum density: p(r,t) < er . This takes into account finite size effects 
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and filling factors. Indeed, since the cells have a finite size, they cannot be 
compressed indefinitely. In this generalized Kellcr-Scgcl model, chemotactic 
collapse leads ultimately to the formation of a smooth aggregate instead of 
a Dirac peak in the standard Keller-Segel model. This regularized model 
prevents finite time-blow up and the formation of (unphysical) singularities 
like infinite density profiles and Dirac peaks. Therefore, the Dirac peaks 
(singularities) are replaced by smooth density profiles (aggregates). 

(ii) Models with anomalous diffusion: Chavanis & Sire [32] studied 
a model with a constant mobility and a power law diffusion coefficient 
D(p) = Dp 1 " 1 (with 7 = 1 + 1/n). This lead to a process of anomalous 
diffusion connected to the Tsallis entropy [33]. For < n < n$ = d/(d — 2), 
the system reaches a self-confined distribution similar to a stable polytropc 
(e.g. a classical white dwarf star) in astrophysics. For n > n^, the system 
undergoes chemotactic collapse above a critical mass (the classical chemo- 
tactic collapse related to the standard Keller-Segel model is recovered for 
n — > +oo) [32]. In the pre-collapse regime, the evolution is self-similar and 
leads to a finite time singularity. A Dirac peak is formed in the post-collapse 
regime. For n = n^, the dynamics is peculiar and involves a critical mass 
similar to the Chandrasekhar limiting mass of relativistic white dwarf stars 
in astrophysics [34] . The case of negative index n < is treated in [35] with 
particular emphasis on the index n = — 1 leading to logotropes. 

In the present paper, we discuss a larger class of generalized Keller-Segel 
models and intcrprete these equations in relation with nonlinear mean field 
Fokkcr-Planck equations and generalized thermodynamics. For illustration, 
we present for the first time a model incorpoating both a filling factor and 
some effects of anomalous diffusion. 



2. The generalized Keller-Segel model 
2.1. The dynamical equations 

The general Keller-Segel model [3] describing the chemotaxis of bacterial 
populations consists in two coupled differential equations 

^ = V ■ (D 2 V P ) - V ■ ODxVc) , (1) 

e-^ = -k(c)c + f(c)p + D c Ac, (2) 

that govern the evolution of the density of bacteria p(r, t) and the evolu- 
tion of the secreted chemical c(r, t) . The bacteria diffuse with a diffusion 
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coefficient Z?2 and they also move in a direction of a positive gradient of 
the chemical (chemotactic drift). The coefficient D\ is a measure of the 
strength of the influence of the chemical gradient on the flow of bacteria. 
On the other hand, the chemical is produced by the bacteria with a rate /(c) 
and is degraded with a rate k{c). It also diffuses with a diffusion coefficient 
D c . In the general Kcller-Segcl model, D\ = Di(p,c) and D2 = D2(p,c) 
can both depend on the concentration of the bacteria and of the chemical. 
This takes into account microscopic constraints, like close-packing effects, 
that can hinder the movement of bacteria. If we assume that the quanti- 
ties only depend on the concentration of bacteria and write D2 = Dh(p), 
D\ = xg{p), k(c) — k 2 , /(c) = A and D c = 1, we obtain 

^=V-(Dh(p)Vp- X g(p)Vc), (3) 

e-Q- t =Ac~k 2 c + \p. (4) 

For e = 0, Eq. (4) becomes the screened Poisson equation 

Ac-fc 2 c= -Xp. (5) 

Therefore, we can identify A: -1 as the screening length. If we assume fur- 
thermore that k = 0, we get the Poisson equation 

Ac = -Xp. (6) 

The generalized Keller-Segel (GKS) model (3) can be viewed as a non- 
linear mean-field Fokker-Planck (NFP) equation [5]. Written in the form 
dtp = V • (V{D(p)p) — x(p)pVc), it is associated with a stochastic Ito- 
Langevin equation 

-= X (p)Vc + V2^)R(i), (7) 

with 

^ )= -M D{p) = - [ P h(x)dx, (8) 
P P Jo 

where R(f) is a white noise satisfying (R(i)) = and (Ri(t)Rj(t/)) = 
SijS(t — t 1 ) where i = 1, d label the coordinates of space. The standard 
Keller-Segel model is obtained when the mobility x and the diffusion coef- 
ficient D are constant. This corresponds to h(p) = 1 and g(p) = p. In that 
case, the stochastic process (7) and the Fokker-Planck equation (3) are sim- 
ilar to the ordinary Langevin and Smoluchowski equations describing the 
diffusion of a system of particles in a potential $(r, t) = — c(r, t) that they 
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produce themselves through Eq. (4). For example, when Eq. (4) is approx- 
imated by Eq. (6), the system becomes isomorphic to the Smoluchowski- 
Poisson system describing self-gravitating Brownian particles [21,22]. The 
steady state of the standard Kcllcr-Segel equation is p ~ ei c . This is similar 
to the Boltzmann distribution p ~ e~*/ T of statistical equilibrium provided 
that we introduce an effective temperature T through the Einstein relation 
T = D/x- In the present study, we shall consider more general situations 
and allow the mobility x(p) an( i the diffusion coefficient D(p) to depend 
on the local concentration of particles p(r, t). This is an heuristic approach 
to take into account microscopic constraints that affect the dynamics of 
the particles at small scales and lead to non-Boltzmannian distributions at 
equilibrium. Indeed, it is not surprising that the mobility or the diffusive 
properties of a particle depend on its environement. For example, in a dense 
medium its motion can be hampered by the presence of the other particles 
so that its mobility is reduced. 



2.2. Generalized free energy and H-theorem 

We define the energy by 



For e = 0, this expression reduces to 



E=-^Jpcdv. (10) 
On the other hand, we define the temperature by 

T=-. (11) 

X 

Therefore, the Einstein relation is preserved in the generalized thermody- 
namical framework. We also set (3 = 1/T. We introduce the generalized 
entropic functional 

S = -Jc(p)dr, (12) 
where C(p) is a convex function (C" > 0) defined by 

C"(p) = (13) 

g(p) 

This defines the entropy up to a term of the form AM+B where M = J pdr 
is the mass (which is a conserved quantity). We can adapt the values of the 
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constants A and B in order to obtain convenient expressions of the entropy. 
Finally, we introduce the generalized free energy 

F = E -TS. (14) 

The definition of the free energy (Legendre transform) is preserved in the 
generalized thermodynamical framework. The free energy is the correct 
thermodynamical potential since the system is dissipative. Thus, it must be 
treated within the canonical ensemble [5,36]. 
A straightforward calculation shows that 

F = -T~ [(-&c + k 2 c-\p) 2 dr- f -±—(Dh{p)Vp-xg{p)Vc) 2 dv. 

J J xg{p) 

(15) 

For e = 0, this equation reduces to 

F = - [ -^—ADh{p)Vp - X g(p)Vc) 2 dr. (16) 

J xg(p) 

Therefore, F < (in all the paper, we assume that e, A, x, D, g, h are posi- 
tive quantities). This forms an H theorem in the canonical ensemble [5,36] 
for the nonlinear mean field Fokker- Planck equation (3). This also shows 
that the free energy F[p, c] is the Lyapunov functional of the generalized 
Keller-Segel model (3)-(4). It is sometimes useful to introduce the Massieu 
function 

J = S-/3E, (17) 

which is related to the free energy by J = —j3F. Clearly, we have J > 0. 
We can now consider particular cases: if D = (leading to T = 0), we get 
F = E so that E < 0. If x = (leading to /? = 0), we have J = S so that 
S>0. 

2.3. Stationary solution 

The steady state of Eq. (3) satisfies F = 0. According to Eq. (15), this 
leads to 

Ac-k 2 c=-Xp 7 Dh(p)Vp~xg(p)^c = Q- (18) 
Using Eqs. (11) and (13), the second equation can be rewritten 

C"{p)Vp - f3Vc = 0, (19) 
which can be integrated into 

C'(p)=/3c-a, (20) 
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where a is a constant of integration. Since C is convex, this equation can 
be reversed to give 

p(v)=F(-0c(r) + a), (21) 

where F(x) = (C')~ l {— x) is a monotonically decreasing function. Thus, in 
the steady state, the density is a monotonically increasing function p = p(c) 
of the concentration. We have the identity 

✓«■) = J^- <*0 

Substituting Eq. (21) in Eq. (5), valid for a stationary state, we obtain a 
mean-field equation of the form 

-Ac + fc 2 c = XF(-f3c + a). (23) 

The constant of integration a is determined by the total mass M (which is 
a conserved quantity). Finally, we note that the generalized entropy (12) is 
related to the distribution (21) by: 

C(p) = - f F- 1 (x)dx. (24) 



Equation (21) determines the distribution p(r) from the entropy S and Eq. 
(24) determines the entropy from the density. 

2.4. Minimum of free energy 

The critical points of free energy at fixed mass are determined by the vari- 
ational problem 

SF + TaSM = 0, (25) 
where a is a Lagrange multiplier. We can easily establish that 

5E = ~ j (Ac- k 2 c + Xp)Scdr - J cSpdr, (26) 



SS = - J C'{p)8pdv. (27) 

The variational problem (25) then leads to 

Ac-k 2 c=-\p, C'(p)=f3c-a. (28) 

Comparing with Eq. (20), we find that a stationary solution of Eq. (3) is a 
critical point of F at fixed mass. On the other hand, we have established 
that 

F < 0, F = ^ dtp = 0. (29) 
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According to Lyapunov's direct method [4], this implies that p(r) is linearly 
dynamically stable with respect to the NFP equation (3)- (4) iff it is a (local) 
minimum of F at fixed mass. Maxima or saddle points of F are dynamically 
unstable. In conclusion, a steady solution of the GKS modcl/NFP equation 
(3)- (4) is linearly dynamically stable iff it satisfies (at least locally) the 
minimization problem: 

min {F[p,c} | M[p] = M}. (30) 

p,c 

In this sense, dynamical and generalized thcrmodynamical stability in the 
canonical ensemble coincide. Furthermore, if F is bounded from below a , 
we can conclude from Lyapunov's theory that the system will converge to 
a stable steady state of the GKS model for t — > +oo. 
Finally, we note that the GKS model can be written 



dt 



, \ t— t &F 

X9KPW-T- 
Sp 



(31) 



where S/dp is the functional derivative. This shows that the diffusion current 
J = —xg(p)V(8F/5p) is proportional to the gradient of a quantity 5F/Sp 
that is uniform at equilibrium ((SF/Sp) eq = —Ta according to Eq. (25)). 
This corresponds to the linear thermodynamics of Onsager. The same result 
can also be obtained from a generalized Maximum Free Energy Dissipation 
(MFED) principle which is the variational formulation of Onsager's linear 
thermodynamics [5], 

2.5. Particular cases 

If we take h(p) = 1 and g(p) = 1/C"(p), the NFP equation (3) becomes 

In that case, we have a constant diffusion D(p) = D and a density depen- 
dent mobility x(p) = x/{pC"{p))- If we take g{p) = p and h(p) = pC"(p), 
the NFP equation (3) becomes 

^=V-(DpC"(p)Vp- X pVc). (33) 



a We note that for the standard Keller-Segel model, or for the Smoluchowski-Poisson 
system, the free energy is not bounded from below. In that case, the system can either 
relax towards a local minimum of F at fixed mass (when it exists) or collapse to a Dirac 
peak [24], leading to a divergence of the free energy F(t) — > -co. The selection depends 
on a complicated basin of attraction. The same situation (basin of attraction) happens 
when there exists several minima of free energy at fixed mass. 
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In that case, we have a constant mobility x{p) = X an d a density dependent 
diffusion D(p) = Dp[C(p)/p]' . Note that the condition D(p) > requires 
that [C(p)/p]' > 0. This gives a constraint on the possible forms of C(p). 

Finally, if we multiply the diffusion term and the drift term in the NFP 
equation (3) by the same positive function A(r, t) (which can be for example 
a function of p(r, t)), we obtain a NFP equation having the same free energy 
(i.e. satisfying an //-theorem F < 0) and the same equilibrium states as the 
original one. Therefore, for a given entropy C{p), we can form an infinite 
class of NFP equations possessing the same general properties [5] . 

2.6. Generalized Smoluchowski equation 

The NFP equation (33) can be written in the form of a generalized Smolu- 
chowski (GS) equation 

^=V.[x(Vp-pVc)], (34) 
with a barotropic equation of state p(p) given by 

p'(p)=TpC"(p). (35) 
Since C is convex, we have p'(p) > 0. Integrating Eq. (35) twice, we get 

TC(p) = P J" ^dp'. (36) 

Therefore, the free energy (14) can be rewritten 

F=±j [(Vc) 2 + fc 2 c 2 ] dv- J pedv + J p J" ^dp'dv. (37) 

With these notations, the //-theorem becomes 

F = -— [ (Ac - k 2 c + Xpfdr - [ —CVp- pVcfdr < 0. (38) 
Ae J J XP 

The stationary solutions of the GS equation (34) satisfy the relation 

Vp - pVc = 0, (39) 

which is similar to a condition of hydrostatic equilibrium. Since p = p(p), 
this relation can be integrated to give p = p(c) through 

(40) 

P 



This is equivalent to 

P'(P) 



77T- ( 41 ) 
p p'{c) 
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This relation can also be obtained from Eqs. (35) and (22). Therefore, 
we recover the fact that, in the steady state, p = p(c) is a monotonically 
increasing function of c. We also note the identity 

p(p) = -D(p)p = Tp 2 

X 

Finally, we note that the relation (40) can also be obtained by cxtremizing 
the free energy (37) at fixed mass writing 5F — a5M = 0. More precisely, we 
have the important result: a steady solution of the generalized Smoluchowski 
equation (34)-(4) *s linearly dynamically stable iff it is a (local) minimum 
of the free energy F[p, c] at fixed mass M[p] — M . This corresponds to the 
minimization problem (30). 

The generalized Smoluchowski equation (34) can also be obtained for- 
mally from the damped Euler equations [5]: 

^+V-(pu)=0, (43) 
dxx 1 

— + (u- V)u = --Vp + Vc-£u. (44) 
at p 

For £ = 0, we recover the usual barotropic Euler equations of hydrody- 
namics. Alternatively, if we consider the strong friction limit £ — > +oo, 
we can formally neglect the inertial term in Eq. (44) and we get £u = 
— iVp + Vc + 0(£ _1 ). Substituting this relation in the continuity equation 
(43), we obtain the generalized Smoluchowski equation (34) with x = 
These hydrodynamic equations (hyperbolic model) have been proposed 
in the context of chemotaxis to describe the organization of endothelial 
cells [6,37-39]. This inertial model takes into account the fact that the cells 
do not respond immediately to the chemotactic drift but that they have 
the tendency to continue in a given direction on their own. Therefore, the 
inertial term models cells directional persistence while the general density 
dependent pressure term — Vp(p) takes into account anomalous diffusion or 
the fact that the cells do not interpenetrate. Finally, the friction force — £u 
measures the response of the system to the chemotactic "force" Vc. Indeed, 
after a relaxation time of the order £ _1 their velocity will be aligned with 
the chemotactic gradient. For £ = 0, Eqs. (43)-(44) lead to the formation of 
filaments that are interpreted as the beginning of a vasculature [6,37-39]. 
These filaments, or networks patterns, are not obtained in the Keller-Segel 
model (parabolic model), corresponding to £ — > +oo, which leads to point- 
wise blow up or round aggregates [7,23] . Note finally, that the GS equation 



C(p) 



= T[C'{p)p-C{p)}. 



(42) 
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(34) can be derived rigorously from kinetic models in a strong friction limit 
£ — > +oo, using a Chapman-Enskog expansion [40] or a method of mo- 
ments [6]. 

2.7. Kinetic derivation of the generalized Keller-Segel 



As discussed previously, the generalized Keller-Segel model (3) can be 
viewed as a nonlinear Fokkcr-Planck equation where the diffusion coefficient 
and the mobility explicitly depend on the local concentration of particles. 
Such generalized Fokkcr-Planck equations can be derived from a kinetic the- 
ory, starting from the master equation, and assuming that the probabilities 
of transition explicitly depend on the occupation number (concentration) 
of the initial and arrival states. Below, we briefly summarize and adapt to 
the present situation the approach developed by Kaniadakis [41] in a more 
general context. 

We introduce a stochastic dynamics by defining the probability of tran- 
sition of a particle from position r to position r'. Following Kaniadakis [41], 
we assume the following form 



Usual stochastic processes correspond to a(p) = p and b(p) — 1: the prob- 
ability of transition is proportional to the density of the initial state and 
independent on the density of the final state. They lead to the ordinary 
Fokker-Planck equation (64) as will be shown below. Here, we assume a 
more general dependence on the occupancy in the initial and arrival states. 
This can account for microscopic constraints like close-packing effects that 
can inhibitate the transition. Quite generally, the evolution of the density 
satisfies the master equation 



Assuming that the evolution is sufficiently slow, and local, such that the 
dynamics only permits values of r' close to r, one can develop the term in 
brackets in Eq. (46) in powers of r — r'. Proceeding along the lines of [41], 
we obtain a Fokkcr-Planck-like equation 



model 



7r(r - r') = u,(r, r - r')a[p(r, t)]b[p(r' ,t)\. 



(45) 




(46) 



dhiK(p) dp 



dp dxj 



(47) 



with 



70) = a{p)b{p) 



K(p) 



a(p) 



(48) 



KpY 
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and 



Ci(r) = - / Viw(r, y)dy, 



Cij (r) = ^ / yiVjw{r,y)dy. 

The moments and Cij are fixed by the Langevin equation 

dv 



— = X Vc + V2DTL(t). 



(49) 



(50) 



(51) 



Assuming isotropy Ci = Jh Cij = D5ij, the kinetic equation (47) becomes 
<9p 



at 



v 



(J + Vj D)7(p)-7(p)^^W 
dp 



(52) 



Now, according to the Langevin equation (51), D is independent on r and 
J = — xVc. Thus, we get 



dt 



7(P) 9p W Vp-X7(p)Vc 



If wc define 



%)=7(P)^^, ff(p)=7&0, 
the foregoing equation can be rewritten 

= V • (DMP)VP - X.9(p)Vc) , 
and it coincides with the GKS model (3). We note that 

hm(p)=C'(p). 

We also have the relations 



Kp) 



h(p) 



VW) C 



-C'(p)/2 



Inversely, 



K(p) 

g(p) = a(p)b(p), C'(p) = In 



a (p) 

Kp) 



(53) 

(54) 

(55) 
(56) 

(57) 
(58) 

(59) 



h(p)=b(p)a'(p)-a(p)b'(p). 



(60) 
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It seems natural to assume that the transition probability is proportional 
to the density of the initial state so that a(p) = p. In that case, we obtain 
an equation of the form 

^ = V ■ (D [b( P ) - pb'(p)] Vp - xpb(p)Vc) . (61) 

Note that the coefficients of diffusion and mobility are not independent 
since they are both expressed in terms of b(p). Choosing b(p) = 1, i.e. a 
probability of transition which does not depend on the population of the 
arrival state, leads to the standard Fokker-Planck equation, or standard 
Kcllcr-Scgcl model (64). If, now, we assume that the transition probability 
is blocked (inhibited) if the concentration of the arrival state is equal to 
an upper bound <7o, then it seems natural to take b(p) = 1 — p/a®. In that 
case, we obtain 

^ = V • (DVp - X P(1 - p/^o)Vc) , (62) 

which will be considered in Sec. 3.5. Inversely, we can wonder what the 
general form of the mobility will be if we assume a normal diffusion h(p) = 1 . 
This leads to b(p) — pb'(p) = 1 which is integrated in b(p) = 1 + Kp where 
if is a constant. Interestingly, we find that this condition selects the class 
of fermions [K = —1) and bosons (K = +1) and intermediate statistics 
(arbitrary K) . The corresponding NFP equation is 

^ = V-(DVp- X p(l + Kp)Vc). (63) 

3. Examples of generalized Keller-Segel models 

In this section, we consider generalized Kcllcr-Scgcl models of chemotaxis 
and show their relation with a formalism of generalized thermodynamics. 

3.1. The standard Keller-Segel model: Boltzmann entropy 

If we take h(p) = 1 and g(p) = p, we get the standard Keller-Segel model 

^=V-(DVp- X pVc). (64) 

It corresponds to an ordinary diffusion D(p) = D and a constant mobil- 
ity x{p) = X- The associated stochastic process is the ordinary Langevin 
equation 

dr 



;t = X Vc + V2DR(t). (65) 
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The entropy is the Boltzmann entropy 

S = - J plnpdr, (66) 

and the stationary solution of Eq. (64) is the Boltzmann distribution 

p = e p c - a -i_ ^ 

The standard Keller-Segel model is isomorphic to the Smoluchowski equa- 
tion with an isothermal equation of state 

p(p) = PT. (68) 



3.2. Generalized Keller-Segel model with power law 
diffusion: Tsallis entropy 

If we take h(p) — qp q ~ l and g{p) = p, we obtain the GKS model 

^ = V • (DVp« - xpVc) . (69) 

It corresponds to a power law diffusion D(p) = Dpi -1 and a constant 
mobility x(p) = X- The associated stochastic process is 

^ = X Vc+y/2D P *?B.(t). (70) 

This model can take into account effects of non-ergodicity and nonexten- 
sivity. It leads to a situation of anomalous diffusion related to the Tsallis 
statistics [33]. For q = 1, we recover the standard Keller-Segel model with a 
constant diffusion coefficient, corresponding to a pure random walk (Brow- 
nian model). In that case, the sizes of the random kicks are uniform and do 
not depend on where the particle happens to be. For q ^ 1, the size of the 
random kicks changes, depending on the distribution of the particles around 
the "test" particle. A particle which is in a region that is highly populated 
[large p(r,t)} will tend to have larger kicks if q > 1 and smaller kicks if 
q < 1. Since the microscopies depends on the actual density in phase space, 
this creates a bias in the ergodic behavior of the system. Then, the dy- 
namics has a fractal or multi-fractal phase space structure. The generalized 
entropy associated to Eq. (69) is the Tsallis entropy 

S = -^-jJ(p«-p)dr, (71) 

and the stationary solution is the Tsallis distribution 

i 

p = 91 [l-(q- l)H3c + a)]^" 1 ' . (72) 



September 29, 2008 15:21 



WSPC - Proceedings Trim Size: 9in x 6in article 



15 

The generalized Keller-Segel model (69) is isomorphic to the generalized 
Smoluchowski equation (34) with an equation of state 

p(p) = Tp q . (73) 
This is similar to a polytropic gas with an equation of state p = Kp 1 (with 
7 = 1 + 1/n) where K = T plays the role of a polytropic temperature and 
q = 7 is the polytropic index. For q = 1, we recover the standard Keller- 
Segel model (64). For q = 2, we have some simplifications. In that case, the 
GKS model (69) becomes 

^ = V ■ {DVp 2 - xpVc) . (74) 

The entropy is the quadratic functional 

S = - J p 2 dr, (75) 

and the stationary solution is 

p=-\{-pc+a), (76) 

corresponding to a linear relation between the density and the concentra- 
tion. In that case, the differential equation (23) determining the steady 
state reduces to the Helmholtz equation. Finally, the pressure is 

P(P) = Tp 2 . (77) 

In the context of generalized thermodynamics, the NFP equation (69) was 
introduced by Plastino & Plastino [42] and the generalized stochastic pro- 
cess (70) was introduced by Borland [43]. When the NFP equation (69) is 
coupled to the Poisson equation (6), we obtain the polytropic Smoluchowski 
Poisson system describing self-gravitating Langevin particles. When the 
NFP equation (69) is coupled to the field Eq. (4), we obtain a generalized 
Keller-Segel model of chemotaxis taking into account anomalous diffusion. 
These models have been introduced and studied by Chavanis & Sire [6,32]. 
For the particular index 713 = d/(d — 2) or q 4 / 3 = 74/3 = 2{d — l)/d, the 
GKS model presents a critical dynamics [34] . 

3.3. Generalized Keller-Segel model with logarithmic 
diffusion: logotropes 

If we take h(p) = l/p and g{p) = p, we obtain a GKS model with a 
logarithmic diffusion 

^ = V-( J DVlnp- XP Vc). (78) 



September 29, 2008 15:21 



WSPC - Proceedings Trim Size: 9in x 6in article 



16 



The generalized entropy associated to Eq. (78) is the log-entropy 

S = J lnpdr, 

and the stationary solution is 

1 



(79) 



(80) 



P = Tr- 
et — pc 



The pressure law is 

p(p) = Thxp. 



(81) 



This is similar to a logotropic equation of state [44] . This is also connected 
to a polytropic equation of state (or Tsallis distribution) with 7 = q = 0. 
Indeed, the logotropic model (78) can be deduced from Eq. (69) by writing 
DWp q = Dqp q ~ 1 \/ p, taking q = and re-defining Dq — > D. When the 
NFP equation (78) is coupled to the Poisson equation (6), we obtain the 
logotropic Smoluchowski-Poisson system. When the NFP equation (78) is 
coupled to the field Eq. (4), we obtain a generalized Kellcr-Scgel model of 
chemotaxis. These models have been introduced and studied by Chavanis 
& Sire [35]. 

3.4. Generalized Keller-Segel models with power law 
diffusion and power law drift: Tsallis entropy 

We introduce here an extension of the GKS model (69). If we take h(p) = 
qpi+v-i anc [ g(pj = pt 1 + 1 ^ wc obtain 



This corresponds to a power law diffusion D(p) = ^^jP 9+Al 1 arj d a power 
law mobility x(p) — XP^ ■ The associated stochastic process is 



Since p M can be put in factor of the diffusion current in Eq. (82), this model 
has the same equilibrium states (72) and the same entropy (71) as Eq. (69). 

For a = 0, we recover Eq. (69) with a constant mobility and a power 
law diffusion. For (jx, q) = (0,0), we recover the logotropic Smoluchowski 
equation (78) provided that we make the transformation Dq — > D. For 
/it = 1 — q, we have a normal diffusion and a power law mobility 




(82) 




(83) 



dp 
db 



V • (DqWp - X P 




(84) 
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For q = 2, we get 

^ = V-(2£Vp-xVc), ( 85 ) 

which has the same entropy and the same equilibrium states as Eq. (74). 
Finally, for q = (making the transformation qD — > D), we obtain 

^ = V ■ {DVp - X p 2 Vc) , (86) 

which has the same entropy and the same equilibrium states as Eq. (78). 
When the NFP equation (82) is coupled to the field equation (4), we ob- 
tain a generalized Keller-Segel model of chemotaxis taking into account 
anomalous diffusion and anomalous mobility. 



3.5. Generalized Keller-Segel models with a filling factor: 
Fermi-Dirac entropy 

If we take h(p) = 1 and g(p) = p(l — p/ctq), we obtain a GKS model of the 
form 

^ = V-(DVp-xp(l-p/*o)Vc). (87) 

This corresponds to a normal diffusion D(p) = D and a mobility x(p) = 
x(l — p/cro) vanishing linearly when the density reaches the maximum value 
= cr - The associated stochastic process is 

dr 



, t = X (l-p/ao)Vc + V2Z?R(t). 

The generalized entropy associated with Eq. (87) is a Fcrmi-Dirac-like en- 
tropy in physical space 

S=-* [{^^ + (l-P-)ln(l-JL)\ drt (89 ) 

and the stationary solution is a Fermi-Dirac-like distribution in physical 
space 

P = ^V^- (90) 

From Eq. (90), we see that, in the stationary state, p < ao- This bound is 
similar to the Pauli exclusion principle in quantum mechanics. In fact, we 
can show that p(r, t) remains bounded by <jq during the whole evolution. 
For <7o — > +oo, we recover the standard KS model (64). 
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An alternative GKS model, with the same entropy and the same equi- 
librium states, is obtained by taking h(p) = 1/(1 — p/cfq) and g(p) = p. 
This leads to 

% = V • (— -Drjo V ln(l - p/a ) - xpVc) . (91) 
at 

This corresponds to a nonlinear diffusion with D(p) = — aa(D/p) ln(l — 
p/cra) and a constant mobility \(p) = X- Equation (91) can be put in the 
form of a generalized Smoluchowski equation (34) with a pressure law 

p(p) = -Ta Hl- p/a ). (92) 

For p <C (Jo , we recover the "isothermal" equation of state p = pT leading 
to the standard Keller-Segel model (64). However, for higher densities, the 
equation of state is modified and the pressure diverges when p — > <t - This 
prevents the density from exceeding the maximum value <jq. 

The NFP equation (87) has been introduced by Kaniadakis & Quarati 
[45] to describe fermionic systems and by Robert & Sommeria [46] in the 
statistical mechanics of two-dimensional turbulence (see also [47]). In the 
context of chemotaxis, the model (87) has been introduced by Hillen & 
Painter [29] and, independently, by Chavanis [5,31]. It provides aregulariza- 
tion of the standard Keller-Segel model preventing overcrowding, blow-up 
and unphysical singularities. The filling factor (1 — p/ctq) takes into account 
the fact that the particles cannot interpenetrate because of their finite size 
a. Therefore, the maximum allowable density is cto ~ l/a d . It is achieved 
when all the cells are packed together. In the model (87), it is assumed 
that the mobility vanishes when the density reaches the close packing value 
(p — > do) while the diffusion is not affected. The alternative model (91) has 
been introduced by Chavanis [5,31]. In that case, the mobility is assumed 
to be constant and the rcgularization preventing overcrowding is taken into 
account in the pressure law (92). We can also multiply the diffusion term 
and the mobility term in the NFP equation (3) by the same positive func- 
tion A(r, t) in order to obtain a more general model with the same entropy 
and the same equilibrium states in which both diffusion and mobility are 
affected by overcrowding. 

3.6. Generalized Keller-Segel models incorporating 
anomalous diffusion and filling factor 

The previous models focus individually on two important effects: anoma- 
lous diffusion (see Sees. 3.2-3.4) and exclusion constraints when the density 
becomes too large (see Sec. 3.5). Here we introduce a mixed model which 
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combines these two effects in a single equation. If we take h(p) = qp q+ ^ 1 
and g(p) = p^ l+1 (l — p/cro), we obtain 

^ = V • (Dqp^-^p - XP » +1 {1 - p/a )Vc) . (93) 

This corresponds to a power law diffusion such that D(p) = 
[Dq/{q + p)]p q+ ^ 1 and a mobility x(p) = XP^(^ ~ p/ CT o)- The associated 
stochastic process is 



§ = X p^(l p/* )Vc + J^Lp^m). (94) 



The generalized entropy corresponding to Eq. (93) is obtained by integrat 
ing twice the relation 

,9-2 

P/<?0 

A first integration gives 



C"{ P ) = T~~^~ 7 ■ (95) 

1 - p On 



C'(p) = qol 1 ®^ ( 



(96) 



where 

rt 



*m(t) = I i dx. (97) 

1 — X 



Therefore, the generalized entropy can be expressed as 

[■p/<*a 

C{p) = qal <f q -2(t)dt. (98) 

Jo 

On the other hand, the equilibrium density is given by p — OQ&~^ 2 [(f3c ~ 
a)/q(TQ -1 ]. Note that these results not depend on p since the term p^ can 
be put in factor of the diffusion current in Eq. (93). 

Let us consider some particular cases, (i) For q = 1, Eq. (93) has the 
same entropy and the same equilibrium states as Eq. (87). (ii) For oo — ► 
+oo, we recover Eq. (82). (iii) For p = and q = 2, we have 

^ = V ■ {DVp 2 - X p(l - p/a )Vc) . (99) 
The generalized entropy is 

S = -2ol j f 1 ~ m f 1 - rfr < ( 10 °) 
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and the stationary solution is 



^0 



1 



o (-0c+a)/2a o 



(101) 



For Co — > +oo, we recover Eq. (76). Dividing the diffusion and the drift 
term by 1 — p/cr 0y we can also consider the alternative model 



dt 



2pD 
1- p/a 



(102) 



which has the same entropy and the same equilibrium states as Eq. (99). 
The pressure law is 



p(p) = -2TO-Q [ln(l - p/a ) - p/vq] 



(103) 



(iv) For (/!, q) = (0, 0) and performing the transformation qD — ► D, or 
directly taking h(p) = 1/p and g(p) = p{\ — p/(Jo), we obtain 



^ = V • (DV In p - X p(l - p/ao) Vc) 



(104) 



This corresponds to a logarithmic diffusion and a modified mobility taking 
into account an exclusion principle through the filling factor. The general- 
ized entropy is obtained from the relation 



C"(p) 



1 



p 2 (i-p/<?vY 



leading to 



C'{p) 



1 

^0 



In 



P 



0£ 



and finally to the explicit expression 



S = - 



1 - — ) In ( — - 1 | dr 



We can consider the alternative model 



dp 
dt 



D 



p(l - p/a ) 



(105) 



(106) 



(107) 



(108) 



with the same entropy and the same equilibrium states. The associated 
pressure law is 



p(p) = -T]n(^-l 
P 



(109) 
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4. Conclusion 

In this paper, we have discussed a generalized class of Keller-Segel models 
describing the chemotaxis of biological populations. We have shown their 
analogy with nonlinear mean field Fokker-Planck equations and generalized 
thermodynamics. We have given explicit examples corresponding to differ- 
ent entropy functionals. In particular, we have considered the case where 
the particles (cells) experience anomalous diffusion and the case where they 
experience an exclusion constraint (volume filling). We have introduced a 
mixed model taking into account these two effects in a single equation (93). 
Of course, we can construct other types of Keller-Segel models which may 
also be of interest. The general study of these models, which combine both 
long-range interactions and generalized thermodynamics, is very rich and 
can lead to a wide diversity of phase transitions and blow up phenomena. 
These nonlinear meanfield Fokker-Planck equations are therefore of consid- 
erable theoretical interest [5] . 
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